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We present a comprehensive study of the slip and flow of concentrated colloidal 
suspensions using cone-plate rheometry and simultaneous confocal imaging. In the 



colloidal glass regime, for smooth, non-stick walls, the solid nature of the suspension 
causes a transition in the rheology from Herschel-Bulkley (HB) bulk flow behavior 



at large stress to a Bingham-like slip behavior at low stress, which is suppressed for 



sufficient colloid-wall attraction or colloid-scale wall roughness. Visualization shows 

ID '. 

q . how the slip-shear transition depends on gap size and the boundary conditions at 

m ■ 

Tj- 1 both walls and that partial slip persist well above the yield stress. A phenomenolog- 



ical model, incorporating the Bingham slip law and HB bulk flow, fully accounts for 
the behavior. Microscopically, the Bingham law is related to a thin (sub-colloidal) 
lubrication layer at the wall, giving rise to a characteristic dependence of slip param- 
eters on particle size and concentration. We relate this to the suspension's osmotic 
pressure and yield stress and also analyze the influence of van der Waals interac- 
tion. For the largest concentrations, we observe non-uniform flow around the yield 
stress, in line with recent work on bulk shear-banding of concentrated pastes. We 
also describe residual slip in concentrated liquid suspensions, where the vanishing 
yield stress causes coexistence of (weak) slip and bulk shear flow for all measured 
rates. 

PACS numbers: 83.50.Ax, 83.60.-a, 83.80.Hj, 83.85.Ei 
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I. INTRODUCTION 



Wall slip is a widespread phenomenon in the flow of various liquids. In Newtonian liquids, 
interest in slip has revived due to its relevan ce for flow in nano-pqrous m edia, microfluidic 



device s and along superhydrophobic su rfaces 



feOOll); 



Bocquet and Barrat 



(J2007I); 



Barrat and Bocauetl (119991 ); 



Neto et al 



Zhu and Granick 



(|2004h ]. The so-called slip-length l s , the 



distance to the wall at which the velocity profile extrapolates to zero, can reach many 
molecular diameters, depending on wettability, and strongly affects flow when it is com- 
para ble to the system dimens i ons. Polym er slip has also received cons i derab le attention. 



e.g. 



fll966h : 



Hatzikiriakos and Dealv ( 



Leger et al 



(119971 ): 



1991 



1992): 



Brochardt and de Genned (Il992[ ); 



Westover 



Mhetar and Archer! (119981 )]. Here changes in the chain relax- 



ation dynamics near the wall govern slip: typically a transit ion from weak slip at sma ll 
flow rate to strong slip (l s exceeding hundreds of micrometers Mhetar and Archer! (119981 )]) 
for large flow rate is s een, driven by chain disentanglement ne a r the wall. In surfactan t 



solutions, recent work 



Sa 



Lettinga and Mannevilld Q2009); 



mon et al. 



(2003); 



Becu et al. 



Becu et al. 



(j200J); 



Manneville et al. 



(120041 ): 



(120071 ) ] has shown more complex behavior, 



where slip and wall interactions are coupled to the shear-banding in these systems. 

The most prominent examples of slip in industrial and daily app lic ations o c cur i n 
flow of compl ex, multi-phase fluids Yoshimura and Prud'hommd (119881 ): iBarnea (119951 ): 



Larson (1 993)1. Meaningful characterization of the bu 



requires proper insight into boundary effects [Buscalll (120101 )] . Over the last decades, 



k flow properties of these systems 



many stu dies of slip in these systems ha ve appeared, e.g. i n par tic ulate (colloidal) sus- 



pensio ns 
(120041): 



Yil mazer and Kalyonl ( 



Kalvonl (I2005h 



Soltani and Yilmazer (2008) 



Russel and Grant (2000): 



(2008); 



1989); 



Jana et al 



Aral and Kalvonl ( 



1995) 



ICohen et al. 



994); 



Hartman Kok et al 



Pcrscl 



Walls et al. 



Princen 



1985 



Katgert et al. 



Meeker et al. 



(1200% 



o et al (11994T) ]. colloidal gels 



Isa et al. 



3uscall et al 



j|2003j); IVaradan and Solomon! ( 



2003); 



(2007); 



1993) 



Gibaud et al. 



Wassenius and Callaghanl (120051) ] a nd emulsions and foa ms i Bertola et al 



(I2004al ]bh: 



Salmon et al. 



(120031 ): 



Denkov et al. 



(2003); 



(120051 ): 



( 120081 )]. Despite this large body of work, it is challenging to gain microscopic 
insight into the nature of slip and understand its dependence on material composition, wall 
properties, and flow rate. Broadly speaking, slip results from depletion of the dispersed 
phase near a smooth wall, giving a low viscosity, high shear boundary layer which reduces 
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the apparent bulk viscosity Q. However, both the structure and the origin of this "layer" 
vary greatly between different systems and the interplay with (non)linear bulk rheology can 
give rise to slip being pronounced either at large or small flow rate. 

Slip at low stress or applied shear rate can occur in dispersions where caging, aggregation 
or 'jamming' leads to a solid-like microstructure and mechanical behavior. The system then 
exhibits a yield stress, a y , below which the micro-structure remains intact but, depending on 
wall interactions and roughness, apparent flow of the material can still be m easured. This 
was rec e ntly studiec 



( I2004al ): 



Seth et al. 



or co ncentrated emulsions and other soft particle pastes 



Meeker et al. 



( 120081 )] in presence of smooth walls. There, elasto-hydrodynamic lubri- 



cation, associated with the particles' deformability, causes a lubrication layer between the 
compressed packing and the wall with velocity dependent thickness. This gives a nonlinear 
relation between the slip stress a and the slip velocity v s , but this mech a nism o nl y occurs for 

(|2004bJ); 



Meeker et al. 



Seth et al. 



yery w eakly or non-repulsive particle-wall interactions 

( 120081 )]. Note that nonlinear lubrication is of strong interest to applications of many other 
soft solids, e.g. hydrogels, where combined wall-net work repulsion and adsorption deter- 
mines the friction properties Gong and Osada ( 201ol )]. 

For suspensions of 'hard' particles, most studies of slip have been performed for non- 
Brownian systems, i.e. at large Peclet number Pe = 775 >> 100, where 7 is the true 
)ulk flow rate and tr the B r ownian relaxat ion time. There, slip occur s both in solid- 
Yilmazer and Kalyon ( 1989 ): Kalyon ( 2OO5I )] and liquid-like Jana et al. ( 1995 )] suspen- 
sions (albeit without consensus on the quantitative behavior), but littl e attention has been 



Hartman Kok et al. 



given to the effect of Brownian motion on slip in colloidal systems 
( 12004] ) found that Brownian motion prevents depletion at low shear rates in dilute colloids, 
so that slip becomes apparent only for Pe > 10. However, for more concentrated colloids, 
'crowd ing' competes with Br o wnian r elaxation. This c a uses s tructural arrest, a glass tran- 
sition Pusey and van Megenl ( 119861 ): Ivan Megen et al\ (119981 )]. which for hard-sphere (HS) 
colloids occurs at a volume fraction of <p g ~ 0.58. The associated change from liquid to 
solid-like rheology [Petekidis et al. (2004)] can significantly affect the slip response. 

The slip and flo w behavior in t h ese co ncentrated HS colloidal suspensions is the subject 

(120081 )] we have given a short account of some of the 



of this paper. In 



Ballesta et al. 



1 We refer here to 'intrinsic' depletion under uniform stress, as opposed to depletion of dispersed phase 
due to shear migration as may arise under non- uniform stress. 
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results, here we prese nt a more extensive study of the behavior. Using rheo-microscopy 



Besseling et al. 



(1200a )]. we show that for 



> 



'<:)■ 



slip becomes dominant in the rheology 



for smooth, non-stick walls. We address the dependence of slip on colloid-wall (van der 
Waals) interactions and concentration and quantitatively describe the transition from slip to 
yielding for different combinations of confining walls via a phenomenological model. Below 



the yield s tress, slip cau s es ful 



microgels 



Meeker et al 



fl2004ah : 



pl ug flow, qu a 



Seth et al 



itatively similar to jammed emulsions and 



(120081 )]. but quantitatively different; above 



a threshold stress, the slip stress increases linearly with the slip velocity (Bingham slip 
response) due to a sub-colloidal solvent layer with velocity-independent thickness. The 
concentration dependence of the Bingham slip parameters shows a direct relation to the 
osmotic pressure and yield stress of the suspensions, which can be understood on a semi- 
quantitative basis but still lacks a full theo retical descripti o n. W e further show that deep in 
the glass regime, shear localization effects [Besseling et al\ (120101 )] accompany the slip-shear 



transition and discuss residual slip in concen trated liquid suspensions 



< 



a feature not 



uncovered previously [Ballesta et all (120081 )] due to its limited effect on the bulk rheology. 

The paper is structured as follows. After a description of the measurement setup and 
suspensions in Sec. [Til we present rheology and local velocimetry results for various hard 
sphere suspensions, colloid-wall interactions and confining walls conditions in Sec. IHIl In 
Sec. \1V\ we describe the phen omenological model for slip and yielding, generalizing the 



results in 



Ballesta et al. 



(120081 )] to account for geometries with various confining walls, and 
in Sec. |V] we compare the predictions with the experimental results. The physical origin of 
the Bingham slip parameters is analyzed in Sec. IVI| along with the effect of van der Waals 
interactions. Sect. IVHI and IVIIII deal with shear localization for <fi > <p g and residual slip in 
liquid suspensions for < <j> g , respectively, and we conclude in Sect. [IX] 



II. SAMPLES AND METHODS 

A. Colloidal suspensions 

We used polymethylmethacrylate (PMMA) colloids of various sizes (radii a = 138 nm, 
150 nm, 302 nm and a fluorescent batch with a = 652 nm, measured by light scattering, 
polydispersity ~ 15%), sterically stabilized with a poly-12-hydroxystearic acid (PHSA) layer 
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Barretl (119741 )] and dispersed in a refractive index (RI) matching solvent mixture of decalin 
and tetralin (viscosity rj s = 2.2 mPas). A few measurements were also conducted using pure 
decalin as solvent. Non-fluorescent (Rl-matched) samples were seeded with ~ 0.5% of the 
fluorescent particles (labeled with nitrobenzoxadiazole, NBD), which served as tracers during 
confocal imag i ng (F ig. QXb)). In these solvents the colloids interact almost like perfect HSs 



Bryant et al. 



(|2002| )]. In the decalin-tetralin mixture the refractive index of the colloids is 
~ 1.50 — 1.51, somewhat larger than the bulk PMMA value (1.49), due to solvent absorbtion 
of the colloids and slight swelling, estimated to be ~ 10 %. 

Batches of different volume fractions <fi were prepared by diluting samples centrifuged to a 



random close packed sediment, with volume fraction 



'sed 



'rep 



. One traditiona 



determ i ne <t> for HS suspensions e mploy s the crystal- fluid coexistence boundaries 



method to 



f|l995h : 



Pusey and van Me gen 



Segre et al. 



(119861 )]. Because of polydispersi t y our suspensions do not 

(120111 )] . We have chosen 



Poon et al. 



crystallize, and has to be determined differently 
the following method. We first measure the mass density of the solvent p s (densitometer: 
Anton Paar DMA 4500). We then measure the density p rcp of the close packed sediment 
after centrifugation by taking a small sample, diluting it by a factor / = 10 by weight 
in solvent and measuring the density p* of this diluted suspension. The density of the 
sediment is then deduced from p rcp = p*p s /[(l + f)p s — fp*}. Samples at a given (p/4> rcp < 1 



for the main experiments are then prepared by adding a solvent mass M s to the known 
weight M t of the stock sediment (giving 0/0 rcp = [1 + j^^j 2 -]' 1 ) an d homogenizing the 
sample rigorously. The densitometry results are given in table HJ p s slightly increases with 
increasing particle size, probably due to a slight increase in RI for larger particles and thus 
a slight change in the Rl-matching composition of the solvent mixture. We emphasize that, 
without knowledge of the colloid mass density p c , this method does not provide </> rm but 
only 0/0 rcp . When mentioning absolute volume fractions, we have assumed 



rep 

rr.p = 0-67, 



Schaertl and Sillescu 



as fou nd in simulations of spheres with a polydispersity of ~ 15% 
(119941 )] (monodisperse spheres give 4> Tcp = 0.64). The value of (p rcp = 0.67 may be used 
to extract the colloid density via (p rcp = (p rcp — p s )/(p c — p s ). The resulting densities 
p c ~ 1.12 — 1.13 g/ml are reduced compared to the bulk PMMA value of ~ 1.2 g/ml, but 
reasonable match what is expected from the above mentioned ~ 10% swelling of the particles 
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We note that the effect of polydispersity and compaction rate on the value of <p rcp is still 



We do not separately take into account the reduction of p c arising from the PHS-layer. The density of 
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sample name 


asm340 


asm247 


asm209 


asml95 


a (nm) 


138 


150 


300 


652 


p s (g.mr 1 ) 


0.9153 


0.915 


0.92023 


0.92543 


PRCP (g-mr 1 ) 


1.0596 


1.0594 


1.07 


1.0644 



TABLE I. Particle size a, density of the index-matched solvent p s and density of the random close 
packed sediment prcp for the different samples. 



Poon et al. 



under study [Hermes and Dijkstral ( 120101 )] and not without controversy, see also 
( 1201 lh ] . Yet, in relative terms, our approach is quite accurate, the main uncertainty being 
a variation 5<j) rcp /<j) rC p ~ 0.005 found between different centrifugation runs via drying and 
weighting. The value of this uncertainty is close to the maximum variation in the 'jamming' 
raction for our polydispersit y, as found in simulations of HSs for different compression rates 



Hermes and Dijkstral ( 120101 )]. but our S(f) rcp (in two tests we performed, using a single batch) 
was obtained under similar centrifugation conditions. Moreover, the sedimentation Peclet 
number Pe s = Ana A Apg c /?>kBT (with Ap = p c — p s and g c the acceleration), lies in the range 
Pe s ~ 0.005 — 1 from the smallest to the largest particles. Thus, 'jamming' effects (which 
may cause 4> se d < 4> rcp ) are strongly limited, except possibly for the a = 652 nm sample. 
Most important, for the analysis of the </> dependence of rheological and slip parameters 



(Sect. [VTT) . we employ the directly measured reduced volume fraction 



the least from possible ambiguities 



Poon et al. 



(2011)] 



rep > 



which suffers 



B. Measurement system 

We use a stress-controlled rheometer (AR2000, TA Instruments) with a cone-plate ge- 
ometry (radius r c = 20 mm, cone angle 6 = 1° unless mentioned otherwise) and a glass 
slide (radius 25 mm, thickness ~ 180 /jm) as the bottom plate. A solvent trap saturates 
the atmosphere around the sample minimizing evaporation. We image the flow with a con- 
focal scanner (VT-Eye, Visitech International) through the glass plate via a piezo-mounted 
objective (oil immersion, magnification > 60 x ), and optics mounted on an adjustable arm 

swollen colloids is then simply p c — (ppmmas~ 3 ) + [1 — s~ 3 ]p s with s ~ 1.1 the swelling ratio. As the 
swelling is not exactly known, we have used bare particle radii, measured in decalin, throughout the paper. 
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FIG. 1. (a) Cone-plate rheometer with transparent plate and optics connected via an adjustable 
arm to the confocal scanner. The cone angle 6 and radius r c are shown, (b) Confocal image of 
an Rl-matched suspension with the a = 652 nm fluorescent tracers showing as bright spots. Scale 
bar: 10 /im. (c) AFM image of the sintered layer of a = 652 nm colloids on the glass slide. The 
color scale marks height variations of ~ 500 nm. Scale bar: 2.5 fim. Inset: confocal image of a 
coated slide on a larger scale, 50 x 50 /mi. 



Besseling et al\ (120091 )] . We take 2D movies in the flow-vorticity (x-y) plane at a frame 
rate < 90 Hz at various distances r from the cone center and at equally spaced depths 
z < h, with z the velocity-grad ient direction and h th e local gapsize (Fig. Ufa)), from which 



we extract, via image analysis (Besseling et al\ ((200^)], the velocity profiles v(z) at different 



r. The typical error bar in the reduced velocity v(z)/v(h) is < 5%, achieved by imaging 



sufficiently large 



Besseling et al. 



disp 



acements, i.e. over timescales of at least a few times l/7 a (see also 
(120101 )] ) . The observation window is ~ 50 x 50 /im 2 , similar to the gap sizes 
where imaging is performed (r = 2-10 mm, h = 35-175 /im). Over the observation window, 
the variation in h (Sh ~ 0.9 /xm) is negligible, 5h/h 1, and the geometry locally mimics 
parallel plates separated by h = rtan(0). Experiments were performed at controlled applied 
shear rate % (using the rheometer's fast feedback), going from high to low rates, unless 
stated otherwise. Stress controlled measurements gave the same results; in particular, we 
do not find differences between 'static' (measured for increasing stress) and dynamic stress 
thresholds for flow or slip (the latter measured on reducing a or j a ), except for possible 
initial hysteresis just after loading of very concentrated samples (Sect. IVTl) . 



C. Wall properties 



We have used the following preparations of the surfaces of the cone-plate geometry: (i) To 
prevent slip, the cone, the glass slide or both can be made rough on a scale similar to or larger 
than the particle radius by spin-coating a <fi ~ 0.3 suspension of a = 652 nm radius particles 
and sintering the resulting dense disordered colloidal monolayer for one hour at ~ 120° C in 
a vacuum oven. This sintering causes adhesion of the colloids to the glass (or metal cone), 
probably accompanied by local redistribution of the PHS stabilizer at these temperatures. 
Nevertheless, the sintering leaves the corrugated nature of the disordered monolayer intact, 
giving a wall roughness ~ 500 nm, see Fig. []Tc) and the inset. The suppress ion of slip 
associated with this rough coating has been evidenced in Besseling et al\ (120091 )] (Fig. 16 
in that paper) and is also shown in Sect. IVHI (ii) Use of the glass slides without coating 
gives a surface which is very smooth on the colloidal scale (local roughness < 1 nm measured 
by AFM ). The slides were used either untreated or cleaned with ethanol, methanol or a 
Piranha solution (98% H 2 SO4:30%H 2 O2 aqueous solution, ratio 7 : 3 by mass ), but these 
cleaning methods did not give systematic differences in slip behavior. We also performed 
temperature-controlled rheological measurements (without imaging) using the smooth glass 
slides as bottom plate, achieved by anchoring the slides on the rheometer Peltier plate using 
thermal paste. The stainless steel surface of the cone (cleaned with acetone) is also smooth 
on the colloidal scale, but interacts differently with the colloids (Sect. IVI Bp . 



III. MAIN EXPERIMENTAL RESULTS 



Rheological studies of concentrated HS suspensions ha v e show n the emergence of a finite 



dynamic yield stress for <fi > 0.57 — 0.58 



Petekidis et al. 



(120041 ): IPham et all (120061 )]. For 



< 



the flow curve, i.e. a plot of the measured stress (<r m ) versus applied shear rate (j a ), 



exhibits Newtonian behavior at low 7 a followed by strong shear thinning at higher j a , while 



for 



> 



suspensions are glassy and have a yield stress a y below which no flow occurs. In 



this case the flow curves exhibit a Hers chel-Bulkley behay ior described by a. 



a y + 0*7™. 



This behavior is geometry independent Pham et al\ (120081 )]. at least for <j) < 0.6, indicating 
that % and the bulk flow rate of the material, 7, are the same. For lar ger <t>, different 
geometries may cause small differences for a ~ a y due to shear-banding jBesseling et al. 
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10 2 10 1 10° io 1 



FIG. 2. Measured shear stress a m versus applied shear rate j a for a = 300 nm colloids at <fi ~ 0.64 
in Rl-matching solvent (□) and decalin (•), using a smooth glass slide and coated cone. Full line: 
fit of the low 7 a branch to Eq. (TjQ), giving rj e e = 38.3 Pa-s, a s = 2.94 Pa. The data for decalin are 
multiplied by a factor 1.5 for comparison, the difference with the Rl-matching solvent is due to 
slightly different </>. 



(l2010l )]. but the bulk cr(7) behavior is nevertheless consistent with a HB form. 

We first illustrate how the rheology for <ft > <p g is affected when the colloid-wall interaction 
is changed. Figure [2] shows the stress measured with the rheometer, a m , versus the applied 
shear rate j a for a suspension in pure decalin as well as in the Rl-matching decalin-tetralin 
mixture, at similar volume fractions. For the Rl-mismatched sample, the flow curve clearly 
shows a yield stress and an overall Herschel-Bulkley response. However, in the Rl-matched 
sample, the flow curve at low j a exhibits a branch with Bingham-like behavior: 

Cm = cr s + 1] cS j a . (1) 

Here a s is a threshold stress below which the (apparent) flow completely stops and rj e e is an 
effective viscosity characterizing the stress increase in this branch. For large 7 a the stress 
attains the same nonlinear behavior as the sample in decalin. As shown below, the small 
7 a behavior marks full slip along the glass slide (due to the strongly suppressed colloid-wall 
vdW attraction) with a vanishing shear rate 7 in the bulk of the sample. 

We next discuss results for different combinations of smooth and colloid-coated surfaces. 
Figure [3] shows flow curves for Rl-matched samples with both walls or one wall coated, 



coated cone, uncoated plate unc oated cone, uncoated plate 




10 2 10 1 10° 10 1 10 2 10 1 10° 10 1 




0.5 1 0.2 0.4 0.6 0.8 1 

z/h z/h 



FIG. 3. (a,b) Measured stress versus 7 a for coated cone and plate (o) and with (a) uncoated 
plate (□) and (ft = 0.59 and a = 138 nm, and (b) uncoated cone and uncoated plate (□) with 
(j) = 0.59 and a = 150 nm. Regime I in (a) and I a in (b) represent full slip along one boundary; 
regime lb in (b) represents full slip along two boundaries; and regime II in (a),(b) mark slip plus 
bulk flow. In (a),(b) the dash-dot curves are Hershel-Bulkley fits with n = 0.5, Eq. ([3]), giving 
a = 6.1 Pa-s 1 / 2 , a y = 5.5 Pa in (a) and a = 8 Pa-s 1 / 2 , a y = 14 Pa in (b); dashed curves are fits 
to the Bingham form Eq. [TJ In regime I in (a), where Eq. (|10p applies, this gives /3 = 8.2 • 10 Pa 
s/m, <j s = 1.8 Pa. In regime I a or lb in (b), where Eq. (|10[) and Eq. (jlip apply, the parameters 
are f3 = 8.7.10 4 Pa-s-m" 1 , o\ = 5.2 Pa, ai = 3.4 Pa; full lines in regime II are the global flow 
curves from Eq. IA2I in Appendix I A II using the above parameters, (c) Normalized velocity profiles 
v(z)/v cone for the suspension in (a) at r = 3 mm with coated surfaces (filled symbols) and at 
r = 2.5 mm with uncoated plate (open symbols) for various j a . (d) v(z)/v cone for the data in (b) 
for uncoated cone and plate, at r = 5.5 mm and various j a . Full lines: linear fits. Dotted lines: 
behavior without slip. 
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Fig.[3]^a), and with two uncoated walls (smooth glass plate and stainless steel cone), Fig.[3]^b). 
With both surfaces coated (Fig. [3]^a), open circles), a Herschel-Bulkley behavior is recov- 
ered, confirming that co ating prevents slip, a s we also directly observed via microscopy 



Besseling et al 



(120071 ); iBallesta et all (120081 )]. If only one surface is coated (squares in 



Fig. E(a)), we observe a similar slip response as described above: below a critical applied 
shear rate 7 QiC bulk shear vanishes in the entire gap (see below and the appendix) and we 
have a Bingham response; we define this as regime I. For j a > 7 QjC the flow curve deviates 
from this Bingham regime and approaches the flow curve obtained with coated geometry, 
this regime is noted regime II. If both walls are uncoated, squares in Fig. E]^b), the flow curve 
exhibits two successive Bingham regimes at low shear rate: I a marking slip at the plate, lb 
marking slip at both surfaces (see Sec. IIVD : the second Bingham regime, lb, has a higher slip 
stress and an effective viscosity r] c s half the value of that in regime I a . Eventually, at high 
shear rate, the curve again tends toward the HB flow curve. 

We now turn to the velocimetry results obtained by simultaneous confocal imaging of the 
flow, Fig. [3^c,d). In Fig. [3](c) the filled symbols show flow profiles with both walls coated. 
These profiles are nearly linear and v (z) reaches zero at the glass plate and the applied 
velocity at the cone, showing that the coating provides a no-slip condition. The open 
symbols in Fig. |3](c),(d) show v(z) for one uncoated surface (c) or two uncoated surfaces (d). 
In both cases v(z) is essentially linear but exhibits slip at the non-coated walls. The profiles 
can be fitted by v = v s + jz, with v s the slip velocity at the plate and 7 the bulk shear rate. 
At the smallest j a , 7 is zero; the suspension sticks to the coated cone and rotates as a solid 
body, slipping over the smooth glass. This causes the apparent flow below the yield stress 
in regime I in Fig. [3]^a) or I a in (b). Figure [3]Jc) shows that plug flow can persist in regime II 
for r = 2.5 mm and 7 a = 0.5 s" 1 , as explained later. With two uncoated surfaces, Fig. EJ^d), 
the velocity profiles for intermediate 7 a show slip at both walls, with solid body rotation at 
a fraction of the cone velocity. This occurs in regime lb and, as shown later, also in regime 
II at small r. Eventually, for applied rates % > 5 s _1 , both in Fig. [3](c) and (d) the bulk of 
the sample starts to yield, 7 > 0. The difference between 7 and 7 a decreases on increasing 
7 a and the flow curve approaches bulk HB behavior. 

The typical behavior described above, with the transition from a Bingham to a HB 
flow regime, is representative for all Rl-matched suspensions with > 0.57 and at least one 
smooth wall. We illustrate this in Fig. [Hand Fig. [51 The former shows the two-step Bingham 
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00.59 

10" 2 10"! 10° 10 1 0-58 0.6 0.62 

y(l/s) <|> 

FIG. 4. (a) Flow curves for two uncoated surfaces, a = 150 nm, and various 4>. The lower branches 
of the full lines are fits to Eqs. (|10|lip . the upper branch is Eq. IA2I for slip at two surfaces. The 
parameters are a\% = [7,5.7] Pa, j3 = 19.0 • 10 4 Pa-s-m -1 , a y = 23.5 Pa, a = 18.0 Pa-s 1 / 2 for 
(j) = 0.61 and a x>2 = [28,16.6] Pa, /3 = 26.3 • 10 4 Pa-s-m" 1 , a y = 82 Pa, a = 50.6 Pa-s 1 / 2 for 
(j) = 0.63 (the 4> = 0.59 parameters are as in Fig. [3b). (b) Ratio of the slip threshold stresses at 
the bottom and top surface, o^/ci, versus 4>, the dotted line is the average value. 

behavior for two uncoated walls, as in Fig. G^b) but now including two other concentrations. 
Both the slip threshold stress (ui for the cone, a 2 for the plate) and a y increase strongly with 
cf). Figure. St^b) shows that the ratio between the slip threshold stress for the two surfaces 
is essentially independent of <fi. In Fig. [5] we show (a — cr s )/Vea versus applied shear rate 
for suspensions with a = 138 nm at various <ft, using a coated cone and smooth glass plate. 
Plotted this way, all curves at small 7 a (regime I) are linear and overlap, following Eq. ([[]). 
In addition to the threshold stress, also the fitted values of the Bingham viscosity show a 
strong ^-dependence, both of which are discussed in Sec. ED 

Figure [6]^a) shows velocity profiles from the experiment in Fig. [3]^a) (<f) = 0.59, coated 
cone, uncoated plate), but at a constant applied shear rate 7 a = 1.1 s -1 , at various distances 
r. For this applied rate, a m is close to its value measured in the absence of slip (see Fig. Eta)), 
but slip is nevertheless visible. For large r, slip and bulk flow are simultaneously present, the 
mechanism governing slip in regime I thus persists when bulk shear is present. The slip gets 
more pronounced for smaller r, which is shown directly in Fig. E^b) via the r-dependence 
of the local shear rate j(r) extracted from Fig. (H^a): 7(7*) decreases strongly with r and 
eventually vanishes completely for r < r y ~ 2.5 mm, with r y a "yielding radius" inside of 
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FIG. 5. Shear stress minus slip stress cr s ((ft) divided by the Bingham viscosity r? e ff(0) for a = 138 nm 
and 4> = 0.585-0.65. a s and rj e g were extracted from fits of the small j a behavior to Eq. ([1]). Dashed 
line: o m - a s = 7 a ry efr . 
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FIG. 6. (a) Velocity profiles for fixed shear rate j a = 1.1 s _1 , <f) = 0.59 and a = 138 nm with 
uncoated glass and coated cone, at selected distances r from the center, (b) Corresponding local 
shear rate versus r. Full line: Eq. ([5]) with h = r tan(0), using Eq. ([6]) and rheological parameters as 
in Fig. E)^a). Inset: corresponding r-dependence of the normalized slip velocity v s = V2 / ('jr tan 9) , 
along with the prediction using Eqs. (JH),© and © with v\ = (full line). 



which no shear is present. These data demonstrate directly that, in presence of slip, the 
stress in a cone-plate, while essentially uniform across the gap, is radially non-uniform: for 
small r, a(r) < cr y , while for larger r, a{r) > o y . 
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The velocity profiles at small % in Figs. [3] suggest that in regime I the suspension moves 
as a solid body, without bulk shear, down to the glass plate. To get direct evidence fo r 

(J2004J); 



this and exc l ude lo cally non-uniform flow as seen in colloidal crystals 



Perks et al. 



Cohen et al. 



( 120061 )]. we imaged the first layers of a suspension of fully fluore scent particles 



Ballesta et al. 



(a = 6 52 nm) in regime I in the Rl-matching solvent. As reported in Fig. 2a of 
(120081 )]. in regime I below the critical applied shear rate 7 ftiC mentioned earlier, we observed 
plug flow with v s = v cone down to the first layer of particles. There is thus no highly sheared 
layer of colloids and the Bingham slip response, Eq. [TJ reflects the behavior of a lubrication 
layer between the first particles and the wall. Due to the smooth wall, the colloids do exhibit 
some layeri ng at the surface, bu t this vanishes completely beyond ~ 4 particle diameters 



(Fig. 2a of 



Ballesta et al. 



(120081 )]). where the structure is that of a fully disordered colloidal 



glass. Even near the wall, no significant ordering was observed within the imaging plane. 



IV. MODEL FOR THE SLIP AND YIELD BEHAVIOR 



The local and global rheology i n the previous section can be describ ed via a model based 



Russel and Grant 



on th e same assumpti on used in [Yoshimura and Prud'hommd (119881 ); 
(120001 ); iKalyonl (120051 )]. i.e. that the stress at the wall due to slip matches the stress in the 
bulk. As suggested by the phenomenology of the smooth wall flow curves, Fig. [3](a),(b) and 
the Bingham form Eq. [1] for the slip branch, the local stress a s n p and the slip velocity Vi of 
the colloids along a wall, labeled by i, can be related via: 



Oslip = &i + PiVi. 



(2) 



Here is a hydrodynamic term reflecting the lubrication between the first layer of colloids 
and the wall and the threshold stress <7j is similar to a Coulomb friction term. Unless 
particles strongly interact with one of the smooth surfaces the mean distance of a particle 
to the surface will mostly be determined by the available free space, which only depends on 
(/>; hence we assume that f3\ = 02 = fi(4>). On the other hand we expect Oi to depend on the 
wall interaction which can differ between the surfaces. A more detailed description of these 
terms will be given later. Below we use the convention i — 1 for the top surface, i = 2 for 
the bottom plate and also assume that <j\ > cr 2 . 
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For the bulk flow, we use the HB form, taking for simplicity an exponent n = 0.5: 

oWk = a y + «7 - 5 . (3) 

Although the measured n can vary (0.35 < n < 0.6) for 1 > 4>/4> r c P > 0.85, as long as n < 1 
the results below remain qualitatively correct. The relative contribution of slip and flow is 
then evaluated via CTbuik = Csiip along with the relation between 7 and V12' 

hj + «i + v 2 = h%. (4) 

Due to the dependence on the gap size h, the results are geometry dependent. We illustrate 
this for infinite parallel plates and the cone-plate geometry. 



A. Infinite parallel plates 

This idealized situation serves as the basis for subsequent analysis of the cone-plate 
geometry. Various cases are possible depending on the relative values of a s , a and a y . We 
first consider the case a < a y . In this regime bulk shear is absent, 7 = 0. For ct 2 < o < 07 , a y 
(regime I or I a ) slip is localized at one surface, leading to a = cr 2 + fih^ a . For a\ , a 2 < a < o y 
(regime lb) the sample slips at both surfaces and we have a = (ctj + 02) /2 + /3hj a /2. 

Whenever a 2 < cry, yielding sets in at a slip velocity t> 2 = (a y — cr 2 )/f3. This defines the 
transition from regime I or lb to regime II, where bulk flow and slip coexists. The relation 
7(70) is obtained by equating Eqs. (j2J) and (J3J), using Eq. fll]). For slip at one wall (v\ = 0, 
o~ 2 < (J y < a < ai), this gives: 



7 = ia ~ % + 7o ^1 - y 1 + (7a - iyU , (5) 

a • Act . a 2 . 

Act = <jy - ct 2 , ly = — , and 70 = (6) 

For (71,2 < cr y < ct, the sample yields and slips at both surfaces. The relation 7(70) is still 
given by Eq. (JH]), with: 

a ^1 + ^2 . 2Act 2a 2 

Act = ct,-^— ,7 y = ^ 7r ,and7o = 7 ^. (7) 

Equation (JHJ) can be written in dimensionless form, giving the following master curve for 
the local shear rate versus applied rate in regime II: 

f = 1 + n - Vi + 20, (8) 
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with r = 7/70 and Q = (j a — %)/jo, in which the h dependence is absorbed. This form 
describes the vanishing bulk shear rate for 7 a — > 7 y and the approach towards the HB curve 
(7 — > 7 a ) for large rate. This can also be described in terms of the slip length l s , see App. IA 31 
entering regime II, l s decreases and approaches l s ~ 7- 1 / 2 for large 7 a . 

B. Cone and plate 

To analyze the cone-plate geometry, the variation in gap size needs to be accounted for. 
Approximating the geometry at a distance r by parallel plates, the local stress is deduced as 
before, with a{r) = aihj tan(6'),7 a ). The measured stress a m is found by integrating a{r) 
over the entire geometry: 

1 f Tc 

o'm = — n / a(r)2irrdr. (9) 
Jo 

We again first consider the case a(r) < a y , i.e. solid body rotation of the entire sample 
(7 = for all r). In regime I or I a , where slip occurs at the plate (02 < cr < 01, a y ) we have: 

a^ = a 2 + 2/3tan(£)r c7a /3. (10) 

which is the Bingham form Eq. ([T|) with r] c g = 2/3 tan(#)r c /3. For complete slip at two 
surfaces <Ji,a 2 < < cr y , regime lb, the solid body rotation rate of the sample is reduced, 
see Eq. ( 1A5I) in App. IA 21 The stress is obtained from a balance of the total stress on the 
bottom plate and the cone (App. IA 2\f : 

a\ = (<n + a 2 )/2 + /3tan(#)r c7a /3. (11) 

The transition I a — > I b occurs when a]^ = a l £ giving Eq. (IA6I) in App. IA 21 

Due to the non-uniform slip velocity and associated nonuniform stress, the yielding tran- 
sition and flow in regime II differ from that for parallel plates. For slip at the bottom plate 
(o"i = 00, vi = 0), when increasing j a , yielding starts when <7 s i ip = <r y , i.e. at a radius 
r y = (<j y — cr 2 )/(/3tan(6 l )7 a ). The applied rate where regime II starts follows from r y = r c : 

7a, c = (<r y -<r 2 )//3r c tan(0) (12) 

For % > 7 a c , the boundary between solid body rotation (r < r y ) and slip and shear (r > r y ) 
moves inward. The measured stress in regime II follows from Eq. [9] integrated over these two 
regions. The result, Eq. (IA2I) in App. IA II is shown for specific parameters in Fig. [TJ along 
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FIG. 7. (a) Full line: HB flow curve with n = 0.5, a = 6.1 Pa.s ~ a5 , and a y = 5.5 Pa. Other 
curves: calculations for different geometries with slip at the bottom plate using /3 = 82000 Pa.s/m, 
a s = 1.8 Pa: cone-plate with 9 = 1°, r c = 20 mm (dashed line, Eqs. ()10IA2p ; parallel plate with 
h = 2tan(#)r c /3 (dotted) and h = tan(#)r c /3 (dash-dot), using Eqs. ()3I5I6[) . (b) Zoom in on (a). 

with that for parallel plates for two different gaps. The transition I-II for the cone-plate is 
smoother than for parallel plates due to the 'mixed' nature of the transition in the former. 

For slip at two surfaces, the transition to regime II is somewhat more complicated due 
to a small difference in a(r) between bottom plate and cone in regime I. Yet, with the 
approximation oq — 02 <C o y the analysis is essentially the same as for slip at one surface, 
see App. IA2[ the transition to regime II is described by Eq. (|T2|) with the substitution 
a 2 ^( ( j 1 + a 2 )/2and/3^/3/2. 

V. COMPARISON WITH EXPERIMENT 

A. Global rheology 

The uniform lubrication layer and solid body rotation in regime I lead to an effective 
viscosity r] c s ~ /3r c tan(#), which can be verified by comparing slip branches for different 
cones. Figure E](a) shows results for coated cones with different r c and 9 and uncoated glass 
plate for the same sample (a = 138 nm, <p = 0.63). The data for r c = 20 mm, 6 = 2° and 
r c = 10 mm, = 4° indeed superimpose for all j a , as all geometry dependence of a m in 
regime I and II (Eq. (IA2I) ) enters via r c tan(6 l ). When plotted versus the velocity at the edge 
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FIG. 8. (a) Flow curves for an Rl-matched suspension with a = 138 nm, <ft = 0.63, using different 
coated cones and smooth plates, (b) a m versus velocity at the geometry edge; symbols as in (a). 

of the geometry, j a r c tan(9), Fig. MJo), all slip branches superimpose. 

The transition from Bingham slip to HB bulk flow behavior in the global flow curves is 
well described by Eqs. (TTOl [XT]) and ( IA2j) . The parameters a, a y and (3 and 0^2 entering 
Eq. flA2j) (for regime II) follow from fits of the global rheology with and without coating. 
Examples of the predicted full flow curves are shown in Fig. [3]Ja) (b) and Fig. HI In Fig. IHa) 
the measured transition from regime I to II for one smooth wall is well captured. In Fig. EJ^b) 
and HI where both walls are uncoated, the transition from slip at one surface to two surfaces 
and the subsequent yielding are also well described. The small discrepancy for <fi = 0.63 can 
be attributed to both a slight reduction of the HB exponent n < 0.5 for the largest <fi and 
the presence of shear localization near the yield stress in this regime, see Sec. IVIII 

B. Local shear rate 

The model can be checked further by comparing flow profiles to the predictions for local 
shear rate and slip velocities in Sec. IIVI To obtain these quantities, v(z) is fitted to v — 
jz + v s . Figure 0(a),(c) shows the reduced bulk rate 7/70 and slip velocity v s at the plate 
from v(z) taken at r = 2.5 mm (shown in Fig. E^c)) and r = 4 mm. As expected from 
Eq. (jSJ), bulk shear starts at smaller j a for larger r, hence for r < r c the solid body rotation 
extends beyond the transition rate 7 ajC ~ 0.135 s _1 given by Eq. (jT2l) . Similarly Fig. |9^b-d-e) 
presents the local bulk rate and slip velocities i>i(7a), ^2 (7a) f° r uncoated cone and plate 
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FIG. 9. Velocimetry results corresponding to the data in Fig. [3j (a),(b) Measured normalized 
local shear rate versus j a for (a) r = 2.5 mm (□) and r = 4 mm (o), and (b) r = 2.5 mm (□) 
and r = 5.5 mm (o). (c),(d) Normalized slip velocity v s /v cone at the glass plate corresponding to 
data in (a),(b). In (e) the normalized slip velocity at the cone is shown. Data in (a),(c) are for 
a = 138 nm, coated cone and uncoated plate. Data in (b),(d) and (e) are for a = 150 nm, uncoated 
cone and uncoated plate. Full lines in (a),(b) for 7 > are given by Eqs. (|5|6p and Eqs. (|5|7|) . 
respectively, with parameters given in the caption of Fig. [3j In (d),(e), the curves in regime lb are 
given by Eq. (1A5I) ; the transition I a — > lb occurs at 7* given in Eq. (1A6I) . The bulk (slip) velocity 
at the second i>i,2 plateau (corresponding to r < r y in regime II) is given by Eq. (|A7|) : In (c) and 
(d),(e) the curves for the largest 7 a (where r > r y ) follow from those in (a),(b) via Eq. (JH). 
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FIG. 10. (a) Local shear rate 7 extracted from the measured v(z), versus 7 a for different 0, a and 
r, see symbols in (b), for slip at the bottom plate only (a = 138 nm data) and at both surfaces 
(a = 150 nm data). The full lines represent Eq. ([5]) with the rheological parameters entering via 
Eq. ^ or Eq. (FT]) with h = rtan(#). Dotted line: 7 = 7 . (b) Normalized local shear rate T versus 
normalized applied rate O for various r, <j> and two particle sizes. The full line shows Eq. (JS]). 



taken at r = 5.5 mm (from Fig. [3(d)) and at r = 2.5 mm. Both the transition from regime 
I a to lb and that to regime II are well described. Due to the small difference between the 
slip stresses a% and 02 for bottom and top, the extent of regime lb is limited. Note that 
v\ and t>2 are essentially constant between 7 ajC and (2a y — o\ — O2)//3rtan(0) — 4 — 7 s _1 , 
which motivated the simplification of the model for two uncoated plates in App. IA 2[ For 
7a ^ 4 — 7 s^ 1 , the bulk shear rate 7 and v 1>2 also match the predictions. 

The radial variation of 7 over the geometry in regime II, shown in Fig. [6](b), is also well 
explained by the model. The line in Fig. [6(b) shows that both the onset of yielding at 
r y ~ (a y — a s )/{f36 A i a ) = 2.35 mm and the r-dependence of the bulk shear rate directly 
follow from the measured rheological parameters. Finally, in Fig. flOlf a) we show 7 versus 
7 a for different <fi and particle sizes along with the model predictions (see caption). Using 
the normalization described below Eq. [HI T = 7/70 and Q = (% — %)/7o, all data indeed 
collapse on the master curve r = 1 + Q — y/l + 2Q, Fig. [TOT b). Overall, both macroscopic 
and microscopic observations thus validate the phenomenological model. 
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VI. ANALYSIS OF BULK FLOW AND SLIP PARAMETERS 



Before discussing the nature of the slip behavior, we first present the bulk rheological 
parameters, measured in detail for the a = 138 nm and a = 150 nm samples using rough, 
coated walls. Figure [TTT a) shows that a y , plotted as function of the normalized distance 
to random close packing, strongly increases towards rep and is reasonably described by 
a y a 3 /ksT ~ 0.01(1 — 4> / (fi rC p)~ 3 ■ This increase of a y is due to tightening of particle cages 
and strongly increasing entropic barriers. The HB exponent n, Fig. ITTT b). is in the range 
0.4 — 0.55 with a slightly decreasing trend for <\> — > cj) rcp Q Further, the HB parameter a 
also strongly increases with 0, Fig. [TTT c). following approximately a oc (1 — 0/0 rcp ) _2 ' 5 — 
(1 — <f) / 4>rcp) n ~ p ■ This implies the following scaling of the bulk HB rheology: 

a = a y + «7 n = a y {4>)[l + A(l - 0/0 rcp ) n ( 7 r B ) n ] = a,(0)[l + A(jr m {<f>)) n ], (13) 

with r m = t.b(1 — 4>/(f)rcp) a characteristic microscopic (ft dependent (in-cage relaxation) 




FIG. 11. Bulk HB parameters for a = 138 nm and a = 150 nm versus 1 — (0 / '(f> rcp ) . (a) Normalized 
yield stress. Line: a y = ao[l — (cj) / <j) rcp )]~ p with p = 3 and oo = O.OlfcsT/a 3 . (b) the HB epxonent 
n. (c) The HB parameter a. Lines: a = A<tqtS[1 — (4> / (p rcp )] n ~ p with A = 10, n = 0.45, p = 3 and 
tb = 0.0285 s and tb = 0.036 s the Brownian times. 



time scale in the glass, tb the Brownian time and A given by A = a/[a y T^(l 



rap J 



These exponents do not necessarily reflect terminal slopes of log(<7 w ) vs loK(7 a ). Somewhat larger values 



(n < 0.65) have also been observed in other PMMA HS suspensions 



Koumakis et al 



torn ] 
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FIG. 12. (a) Normalized lubrication parameter fiajr\ s versus 1— 4> / 4> rcp . Symbols: data for different 
particle size (see (b)); full line: result based on Eq. (|17p. Dotted line: result based on Eq. ([P. 



and the explanation in the text. Inset: (•) A versus 1 — 4>/4>rcp from a numerical evaluation of 



Eq. (fl6j) . The full line represent A oc il/Ho- (b) Normalized slip stress a s a 3 /ksT versus 1 — Y , ^ r cp, 
the dotted line is a s = 0.45II, the full line represents Eq. (|18p with A = 0.005 and m = 2.5. 



While this phenomenological behav ior has no solid 



and has been used succesfully in 



Besseling et al 



t heoret ical basis, it is physically plausible 
( 20lu[ )]. We do not further pursue its 



interpretation here, but instead turn to the behavior of the slip parameters. 



A. (p dependence for non-stick walls 

As mentioned, for smooth, non-stick walls the Bingham slip parameters /3 and a s in 
Eq. ([2]) show a strong increase with <fi. We show the parameters extracted from the rheology 
in Fig. IT27 a).(b). where in (a) the lubrication parameter has been normalized by the solvent 
viscosity and particle size, while in (b) the stress is normalized by fc^T/a 3 . We observe 
a characteristic increase of the parameters for (j)/d rC p — > for the different particle sizes, 
except for deviations for the a = 652 nm particles l. 

We now attempt to rationalize this behavior in terms of the physical properties of the 
suspension and the (hard) interaction with the wall. We first discuss the lubrication pa- 
rameter /3. The linear increase of stress with velocity, Eq. ([2]), implies that /3 is governed 
by hydrodynamic friction between the first layer of particles and the wall, with an effective 
lubrication layer thickness independent of velocity. For a single particle of radius a, centered 
4 The latter might be partly due to a somewhat reduced sediment volume fraction cf> se d < <j> r cp from 
centrifugation of this batch (as discussed in Sect. Ill Ap . thus slightly overestimating of 
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at a distance a + 5 from the wall, moving with constant velocity v in a solvent of viscosity r] s , 
the drag force it experiences is rj s vaf(5), where f(S/a) ~ 18.1 — 101n(5/a) for S <C a. f(6) 
reflect s the wall induced hydrodynamic reduction of the particle mobility Goldman et al 
( 119671 )] for no-slip boundary conditions. For a distribution of particles with concentration 
n(5) moving with respect to the wall, the mean stress on the wall follows from integrating 
over 6. The lubrication parameter /3 = a/v may thus be written in normalized form as: 



(3a/r] s = aa/r] s v 



a 2 f{S)n(5)d5. 



(14) 



We assume that beyond the first layer the fluid and colloids move together, hence the integral 
can be cut off at 5 = a. With the further assumption that solvent flow in the lubrication 
layer does not change the distribution n(5) and that colloid interactions can be ignored, 
may be evaluated from the "equilibrium" distribution n e (5). It was demonstrated in 



Henderson and van Swoll (119841 )] that the contact value n e (0) follows n e (0) = U/ksT with 



I I the osmotic pressure of the suspension at r e st. While H(4>) is unc e rtain for <ft > 



Hermes and Diikstral ( 



used form is 



201 



0) 



Phan et al. 



Brady and Vivid ( 



199~1 ]: 



(119961 ). iTokuyama and Teradal (120071 )]). a widely 



n = 2.9n /(l - 0/0 rC p), with n = 30fc B T/(4vra 3 ). 



(15) 



Simulations for 6 < 0.5 in 



Henderson and van Swoll (119841 )] further showed that in the first 
layer n e decreases as n(5) = n e (0) exp(— 3X6(1 + 5/a + 5 2 /3a 2 )/a), with A a ^-dependent 
parameter. Employing this form also for the colloidal glass, A(0) can be extracted by 
equating <p in the first layer to the bulk value: 

2vra 2 



n(5)d5 = <p ~ a / n e (0) exp(-X5/a)d5 ~ n e (0)a 6 /X 



(16) 



3 , 

Numerical solution of the l.h.s. of Eq. [161 gives X(<p) oc n(^)/(/c B T/a 3 ) (inset to Fig. IT2T a)). 
as also confirmed by the scaling in the r.h.s. of Eq. [16j Using A(0) we then calculate (3 from 
Eq. [TH The result, shown by the dashed line in Fig. H2T a). is finite for all <ft and diverges as 



'rcpi 



approximately as (3a/r] s oc — ln(l 



rep 



). While this qualitatively accounts for 



the data, and Eq. [T41 matches the experimental oc r/a -1 scaling (Figs. [TBI and [T2Ta)). the 
predicted f3(4>) clearly does not properly describe the experimental results. 

A possible explanation for the discrepancy is that the chosen form for f(S) only applies 
for a single particle. However, while colloid interactions are known to limit the wall induced 
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reduction of diffus i on in concentrated colloidal liquids [Michaifidou et al\ (120091 )]. recent 



simulations 



Swan! (120101 )] show that for small 5 the logarithmic f(5) nevertheless holds. 



An alternative explanation is that the "equilibrium" form n e does not correctly represent 
the near wall particle distribution. This can be due to the fact that already at rest the 



structure is out of equilibrium (since 



> 



possib ly combined with 



change in concentration in the wall layer as observed in Dullens and Kegel 



he lay ering or a 
(120041 )] . Equally 



likely is that n(S) differs from n K due to the actua l flow present in the slip layer, similar 



to dilute systems 



Polverari and van de Yen 



(119951 )]. This non-equilibrium effect can be 



quantified by comparing the shear rate in the lubrication layer with the inverse timescale 
l/r m for "cage" exploration (discussed at the start of this section), via a "wall" Peclet 
number Pe^ = r m v s /^ = r m (a — cr s )/r] s (£ is a mean layer thickness discussed below). Here 
r m is estimated usin g the m e an fre e particle space (also described below) and the short time 



diffusion coefficient 



Bradv! (119961 )] as r m = 3r B [l 



rcp ], in line with r m inferred from 



the 0-scaling of the bulk rheology at the start of this section. Using the Brownian time 
t b > 30 ms and a — a s > 0.2 Pa, we have Pe w = 0(1) and larger for our data in t he slip 
regim e. Thus, during slip n(5) is indeed expected to differ considerably from n e Brady 
(119931 )]. but further measurements are required to confirm this. 

Empirically, the bulk of the experimental data are well described by the form: 



(3a/ Vs = a/Z ~ 0.9/(1 



rcp J 



n/3.2n , 



(17) 



where the l.h.s. defines the mean lubrication layer thickness £ oc 1/11. In fact, the 
mean spacing (s) between colloids in the bulk is estimated as (s) = a((0 rcp / '(ft) 1 / 3 — 1) ~ 
(a/3)(l — (0 / (prep)) , hence the empirical form implies £ ~ 3.4(s). Using n e to calculate £ 
gives similarly £ = a 2 J Q a 5 n(5)d5 ~ 2an /n ~ 2(s), but this form can also result from 
other distributions, different from the 'equilibrium' one. Overall, the analysis suggests that 
the non-equilibrium behavior in the lubrication layer underlies the behavior of j3((p). Note 
that our /3(<p) matches the predicted dependence of the h igh frequency viscosity or inverse 
short time diffusion constant in concentrated suspensions Bradyl (119961 )]. for which, to our 
knowledge, no experimental verification yet exist for <p — > <p rcp . 

Next we discuss the behavior of the slip threshold cr s . As seen in Fig. [T2Tb). a s oc 1/a 3 , 
suggesting a relation to the osmotic pressure. In this context, a s and II may be naively 
linked via a phenomenological Coulomb friction mechanism using a s = fill with \i the 
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FIG. 13. (a) Measured stress a m versus j a for a suspension with a = 138 nm, eft = 0.61 in RI- 
matching solvent, immediately after loading, going from small to large j a (A), followed by large 
to small 7 a (y), again small to large(o) and finally large to small 7 a (o). (b) The power m and (c) 
a s , both obtained from a fit of the low -y a branch to a m — a s oc j™, versus run number. 

friction coefficient. Comparing this with the data, we see that a value \x ~ 0.45 can describe 
<j s at intermediate 1 — 0/0 rC p- Furthermore, recalling the experiments with two uncoated 
surfaces in Fig. HI we found that the ratio between a s at the bottom plate and the cone is 
insensitive to 0, Fig. H^b), which might support the Coulomb friction scenario. However, 
the overall dependence of the data in Fig. IT27b) is inconsistent with a s oc II for constant 
fi. The data are considerably more scattered than for j3(<p) and also exhibit unsystematic 
variation of a s within a factor 2 for different glass cleaning methods. While preventing a 
precise description of the divergence, the data are reasonably described by: 

a s a 3 /k B T ~ A{1 - 0/0 rcp )" m , (18) 

with A ~ 0.005 and m ~ 2.5. This dependence mimics the behavior of the yield stress, 
a y a 3 /ksT ~ 0.01(1 — 4>/4> r c P )~ 3 - Thus, o s ja y = O(0. 1) and decreases weakly with as 
o s joy ~ 0.5(1 — (p/<f)rcp) ' 5 ■ Therefore cx s (0) may represent the stress required for the pre- 
sumed reorganization of the particle distribution when slip sets in and the similarity between 
cr s (4>) and <7 y ((fi) could be connected to a similar change of the particle distribution due to 
'cage' breaking for a = a y . 

A definite interpretation of a s is thus still lacking. We conclude by showing that never- 
theless normal stresses in the system do seem to affect a s . For concentrated samples, the 
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FIG. 14. Mean squared displacement of a = 652 nm particles versus time in dilute suspensions in a 
decalin-tetralin mixture far from the glass ( ■), close to the glass (•) and close to the glass in decalin 
(a). The lines represent diffusive behavior Ax 2 = 2Dt with D = 0.139 /im 2 /s, D = 0.062 /mi 2 /s 
and D = 0.014 /im 2 /s from top to bottom. Inset: trajectories at z ~ a of particles in decalin (the 
stuck particle and short trajectory) and a single long particle trajectory in decalin-tetralin (•). 

rheology of the system may show a transient behavior after initial sample loading, associated 
with local shear thickening behavior during loading. This is shown in Fig. [12] where a sample 
is submitted to repeated low-to-high-to-low shear rate cycles after loading. Fitting the low 
shear part of the first flow curve (A) to a m — cr s oc 7™, we find m ~ 0.65, i.e. in between the 
exponent 1 for Bingham slip and 0.5 for bulk flow, indicating a mixture of shear and slip. 
For the first run a s is significantly larger than for the next cycle(s), where both a s and m 
reach a constant value o s ~ 2 Pa and m ~ 1. Assuming that local shear thickened regions 
lead to relatively large local normal stresses, increasing a s , the reduced a s in repeat runs is 
consistent with flow induced relaxation of these local normal stresses. 



B. Effect of wall interaction 



As discussed, the slip response for HS glasses strongly depends on the colloid-wall vdW 
interaction Q As shown in Fig. [2J for sufficient attraction slip vanishes, implying that the 
5 Note that, generally, attractive colloid-wall forces alone do not guarantee elimination of slip. For e.g. 
flocculated gels, with a heterogeneous network of aggregates, a smoo th, attractive (or even particle coated) 



wall can act as an easy 'fracture' plane resulting effectively in slip 



Buscall et al 



(|19931 )]. In these cases 



larger scale wall roughness - of the order of the aggregate size - is needed to prevent slip. 



0.7 



H 0.6 
m 

CO 

S 0.5 



(a) 



10 20 30 40 50 



12000 
9000 
p^ 3 6000 
3000 




(b) 




10 20 30 40 50 
T(°C) 







TfC) 



10 20 30 40 50 



T( U C) 



FIG. 15. Temperature dependence of (a) the normalized Bingham slip stress and (b) Bingham 
viscosity rj e ff(T) normalized by the temperature dependent solvent viscosity ij s (T). Data are for 
4> = 0.59, a = 138 nm using a coated cone and smooth glass. The inset to (b) shows the un- 
renormalized data r] e g(T). 



slip threshold stress a s > a y . To characterize directly the attraction, we imaged dilute 
suspensions of the a = 652 nm particles and analyzed the near-wall motion both in RI- 
matching and RI- mismatching solvents, without flow. In Fig. [H] we show the mean squared 
displacements (MSD) for the two cases, for particles imaged at the surface (i.e. z ~ a 
within the microscope z— resolution) and t h e MSD away from the su rface. As expected, 



Goldman et al. 



(1967) 



Sharma et al. 



(J2008|); 



Michailidou et al. 



200% } the MSD is smaller 



close to the wall than in the bulk. Moreover, with Rl-matching the near wall MSD is five 
times larger than in the RI- mismatching solvent; in the latter case particles are (temporarily) 
stuck to the surface, evidencing the vdW attractions. With Rl-matching the mobility is 
enhanced, vdW forces are reduced and insufficient to stick particles to the surface. 

Even though the vdW interaction is strongly reduced by Rl-matching, it cannot be com- 
pletely suppressed. To study the role of remaining vdW forces on the slip stress a s (see 
Eq. [1]), we have measured flow curves on the same sample, with smooth glass and coated 
cone, at various temperatures T. Changing T changes the RIs (mainly of the solvent) and 
might thus be observable in the dependence of a s on T. For all temperatures measured 
(5 — 45 C°) the flow curves exhibited the Bingham slip branch (data not shown), from 
which we extracted u s (T) and the effective viscosity fj e s{T) via Eq. [TJ Figure [ToT a) shows 
that the normalized slip stress, a s a 3 /kBT, is nearly temperature independent, but exhibits 
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a small drop of ~ 10 % for T = 25°C. For decalin 8tUr = —4.4 • 10~ 4 K -1 , introducing a 
temperature dependent solvent refractive index ur in the calculations presented in App. Q3] 
can be used to obtain a T-dependent particle-wall interaction. For what interests us here, 
there is a temperature interval where »pmma < "^solvent < ^R.giass and thus colloid-wall 
interactions may become slightly repulsive (see Table 3 below). The observed minimum of 
a s can qualitatively be associated with such slight repulsion, although a calculation suggests 
that this should occur around 40°C rather than 25°C as observed experimentally (figure 
[ToTa)). We attribute the difference to the approximations made in the calculation of the 
vdW interactions (App. [B]). Overall, we conclude that vdW forces do affect a s but, in the 
range of RI- mismatch considered here, have only a modest effect. The data for the Bingham 
viscosity r] eS (T), inset to Fig. [T5Tb). also provide useful information. While r} e g(T) exhibits 
a clear temperature variation, when normalized by the temperature dependent solvent vis- 
cosity, r] s (T) no T dependence is detected (Fig. [ToT b)) suggesting that the viscous slip is due 
to a lubrication layer of pure solvent between the colloids and the wall. 

Rheology experiments along with near-wall motion measurements evidenced that suffi- 
ciently strong vdW attraction between colloids and walls suppresses slip in the HS glasses 
even with smooth walls. Changing the index matching influences both the vdW interaction 
between two particles and the particle-wall interaction. We calculated the interaction V pw 
between the wall and a colloid of radius a separated by a distance 5, and the interaction 
energy V pp between two colloids separated by a distance 25 using the formulas in App. El 
The parameters for each solvent are given in table [Til The resulting particle-particle and 
particle- wall interactions for the various a and the different solvents are shown in table II III 
using 5 = 10 nm as thickness of the steric layer. In decalin, particle-particle vdW attractions 
are weaker than fc^T/10 and thus can be neglected. However, the particle- wall (glass) vdW 
attraction is of the order of ksT for the larger particles explaining the significant tendency 
of these particles to stick to the surface, although such attractions are clearly weaker than 
/cgT/10 for smaller particles. Finally, the vdW attraction between colloids and the metallic 
cone is stronger than k^T for all particle sizes suggesting that a layer of stuck particles 
at the cone should be expected. In decalin-tetralin mixture all attractions are at least an 
order of magnitude smaller than ksT. Hence in index matching solvent particle-particle and 
particle-wall interaction are reduced to their hard sphere counterparts. 

Still, even with a good index matching, some residual van der Waals forces are present 
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Component 


decalin 


decaline-tetraline 


PMMA 


swollen PMMA 


glass 




1.47 


1.51 


1.49 


1.51 


1.523 


e 


2.43 


2.63 


2.6 


2.6 


3.4 



TABLE II. Dielectric permittivity e and index of refraction n B for the solvents, glass, and PMMA. 

(see table IHIl ). Moreover, in decalin-tetralin mixtures, the particle-glass surface interaction is 
positive, which denotes a repulsive force. This is due to the fact that in this case ei < e 3 < e 2 
(see table HT])- Even if such a repulsive interaction should enhance slip, the small energies 
involved allow us to mostly neglect this effect. Two other factors that can affect the slip 
parameter by modifying the friction coefficient are the variations of the index-matching 
solvent composition and the polymer stabilizing layer between different batches and particles, 
and the glass plates that are replaced in each measurement. These factors may introduce 
some experimental uncertainty resulting in a variation of a s by up to a factor of 2 between 
two different experiments. Thus, the friction coefficient may vary both for different particle 
sizes and between different experiments with the same sample. 



sample name 


asm340 


asm247 


asm209 


asml95 


a (nm) 


138 


150 


300 


652 


K c (k B T) 


-0.011 


-0.012 


-0.032 


-0.082 


V p t c (k B T) 


-0.099 


-0.11 


-0.24 


-0.55 


V™ x (k B T) 


-67 e-7 


-76e-7 


-20 e-6 


-51 e-6 


V™ x (k B T) 


0.00097 


0.0011 


0.0024 


0.0054 


V p d l c steel (k B T) 


-4.58 


-5.05 


-11.1 


-25.4 


V™% eel (k B T) 


0.0076 


0.0084 


0.018 


0.042 



TABLE III. van der Waals interactions for different particle sizes a, surfaces and solvents: particle- 
particle interaction in decalin V^f°, particle- wall interaction in decalin V^f, particle-particle in- 
teraction in decalin-tetralin V^ tx , particle- wall interaction in decalin-tetralin V^ lx , as well as the 
particle-steel cone interaction for decalin (V^ c steel ) and for decalin-tetralin (V™ lx teel ). 
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VII. SHEAR LOCALIZATION 



So far we have focused on results for which, at a given r (i.e. gap size) the shear rate 
'j(z) is essentially uniform over the gap when the suspensions start to yield. However, in 
experiments where slip is suppressed by coating both cone and plate, we have observed shear 
banding for small shear rates (near a y ) and large 0. An example is shown in Fig. IToT a). 
For large rates, v(z) is approximately linear, but for 7 < 3 s _1 the profiles become strongly 
nonlinear, with shear localization detected n ear th e walls and a vanishing shear rate in the 



bulk for the smallest j a . In 



Besseling et al. 



(120101 )] . we have shown that this behavior, with 
a continuous variation of j(z) over the gap, can be explained by very small concentration 
gradients (5(f) < 0.003), caused by a d ilation-like flow instability due to shear concentration 



coupling (SCO 



Schmitt et al. 



(11995|)]. This is qualitatively different from other soft glasses 
Besseling et al\ ( l2010l )] and cannot be explained using earlier models for heterogeneous 

Bocquet et all ( 120091 )]. This instability 



glassy flow, e.g. involving specific wall rheology 
and the associated nonlinearity in v(z) sets in below a typical rate 7 C (0) which becomes 
appreciable only for large 0. 



While localization is most easily observed for coated walls, where the average bulk shear 
rate (7) = [v(z g )—v(z = 0)]/z g equals %, we have also detected nonlinear profiles for smooth 
walls near the slip to yield transition, (7) — > 0, for large 0, Fig. [ToT b). For j a = 3 — 5 s _1 
the suspension has started to yield, i.e. (7) > 0, but the j(z) is strongly reduced (enhanced) 
near the smooth (rough) wall compared to (7). Indeed, for 7 a = 5 s" 1 the mean shear rate is 
(7) ~ 0.2% = 1 s _1 , similar to the value of 7 a below which banding becomes significant for 
rough walls, Fig. fT6T a). Due to such nonlinear flow near yielding, it is difficult to determine 
o y very precisely for 0/0 rcp > 0.94. The slightly reduced concentration <p(z) < in the 
fluidized bands allows flow for a (very) slightly below the average o y . Thus the determined 
o y for large <fi may be slightly underestimated, preventing exact calculation of j(r) for 
large where 7 C is large. Therefore, we have only analyzed the slip to yield transition for 
[v(z g ) — v(z = 0)}/z g > 7c, i.e. where the induced bulk flow is essentially linear. 
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FIG. 16. (a) Velocity profiles for = 0.62, r = 5 mm and both surfaces coated, for various j a . (b) 
Same as (a) but for coated cone, smooth plate and r = 5.5 mm. 



VIII. SLIP BELOW THE GLASS TRANSITION 



In our earlier work Ballesta et all (120081 )]. we reported that for concentrations well below 
the glass transition, < 9 , the flow curves and velocity profiles showed no indication of 
slip or shear-banding independent on the surface roughness. Figure [17] illustrates this for 
= 0.52 where the flow curves for coated and smooth surfaces are essentially identical and 
the flow profiles (FigJTTJ (b)) are very close to linear for all j a . However, recent velocimetry 
data for suspensions closer to the glass transition clearly reveal slip at the lowest applied 
shear rate, although the effect is not detectable in the rheology of these samples. 

Figure [18] shows the flow curve and velocity profiles for a = 0.535 suspension with coated 
cone and smooth glass surface. In the flow curve, inset to Fig. H8T a). the low j a Newtonian 
behavior could not be resolved; only the shear thinning behavior could be detected. However, 
the velocity profiles in Fig. [T8T a) clearly show slip, but with a finite bulk shear rate 7 > 
(no plug flow), similar to what is observed in glassy samples for a > a y . We extracted the 
dependence of 7 on v s for various gap sizes. The results in Fig. [T8T b) show that 7 = Cv 2 . 
This is consistent with the model presented earlier, but taking into account the fact that in 
the concentrated liquid regime (here = 0.535) the yield stress is absent. Using a s = 0, 
we obtain a ~ (3v s and a = aj ' 5 giving C = ((3 /a) 2 . From the imaging data in Fig. [181( b) 
we find C = 0.13 • 10 12 s/m 2 at this 0. To compare this with the rheology, a fit of the 
flow curve gives a = 0.177 Pa.s a5 , while follows from extrapolation of the relation /3 = 
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FIG. 17. (a) Flow curve for an Rl-matched suspension with a = 138 nm, eft = 0.52, with smooth 
(□) and rough walls (full line), (b) v(z) for smooth walls at r = 5.5 mm and for various j a . 




FIG. 18. (a) Normalized flow profiles at r = 12.5 mm for <j) = 0.535 at various j a using a smooth 
cone and smooth glass. Inset: measured flow curve, (b) Local shear rate 7 versus slip velocity v s 
at various r. Full line: fit to a power law 7 oc v v s with v = 1.92 ± 0.3. Dotted line: 7 = Cv 2 s with 
C = 0.13 • 10 12 s/m 2 . Magenta discontinued line: 7 = Cv 2 s with C = (/3/a) 2 (see text). 

O.9r] s a _1 (0 rcp — 0) _1 in the glassy state. This yields (3 ~ 5.8 • 10 4 Pa. s/m at this <f>, giving 
((3 /a) 2 = 0.11 • 10 12 s/m 2 , in very reasonable agreement with the data. The model can thus 
also describe residual slip of liquids, assuming a s>y = 0. 

We can compare the slip behavior of glasses and liquids further as follows. The rhe- 
ology of liquid-like suspensions at low 7 is characterized by a Newtonian flow (a = 7707 
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with rjo = 770(0) the zero shear viscosity), which evolves towards nonlinear shear thinning 
behavior at higher shear rates ( a oc 7™). As a result of the low shear rate Newtonian 
behavior, any applied stress results in shear at all r's for any such small j a . At small 7, the 
reduced apparent viscosity is rj m = r}o/(l+rjo//3h) (in plate-plate approximation^). A rough 
estimate of the difference between slip and no slip measurements can be obtained by us- 



ing th e phenomeno 



e.g. 



Meeker et al 



ogical form for the divergence 770(0) on approaching the glass transition 



(11223)]): 770 = 77,(1 - 



of = O^a- 1 ! 



■'rep 



to 



< (j) g . For 



and the earlier mentioned extrapolation 
0.52, a = 138 nm and h = 50 jum this 



leads to i] m = 0.958877 . A ~ 4% difference is roughly within the exp e rimen tal uncertainty, 
which explains why we did not note slip at low in [Ballesta et all (120081 )]. However, as 
approaches <p g , r] diverges while r\ m tends towards (3h and slip becomes apparent. For 
intermediate and large 7, where the rheology can be approximated by a o c 7 n wit h re ~ 0.5 
(as stated above, and in agreement with a semi-empirical expression in iKriegerl (119591 )]). 
a calculation of the flow curves with and without slip shows that in all cases the relative 
difference between r] m and rj increases with and decreases as j a increases, similarly to 
the glassy state. Since in the low-shear Newtonian regime the difference is already small 
(except very close to 9 ), the presence of slip is even more difficult to detect in the nonlinear 
rheology of liquids at larger rate. However, for > 9 the presence of a yield stress and plug 
flow leads to a stress difference ~ a y regardless of the geometry, which explains why slip is 
easily detected in the rheology of HS glasses both for low and intermediate shear rates. 



IX. DISCUSSION AND CONCLUSIONS 



The rheology and velocimetry results and the modeling we presented clearly show that 
local particle-wall interactions, the character of the boundaries, the geometry and de- 
pendent nonlinear bulk rheology strongly affect the measured rheology of concentrated HS 
suspensions. In general for yield stress fluids, if a flow curve exhibits a 'kink' and a stress 
drop at low j a (with a power law 7™) that is absent when rough surfaces are used, one can 
conclude that (i) the sample slips at one or both surfaces and (ii) locally the transmitted 
stress is proportional to v™, independent of the geometry. The detailed phenomenological 

6 In cone-plate, setting h — rtan(#) and integrating as in Eq. [9l the refined expression for the apparent 
viscosity of a Newtonian liquid with slip is rjm/vo — 1 ~ 2x c [l — x c ln(l + a:^ 1 )] with x c = 770/ (j36r c ). 
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model we presented shows that, from a well characterized slip law a(v s ) and bulk rheology, 
the overall flow curve and local flow profiles can be accurately predicted (if more complex 
behavior such as shearbanding and possible non-stationary behavior can be ignored). 

For our HS glasses with smooth non-stick walls, a velocity independent, but ^-dependent, 
lubrication layer forms, leading to m — 1, but vdW attractions, in non-RI- matched suspen- 
sions, easily suppress slip, leading to a slip stress cr s > a y . Standard (non-imaging) rheology 
experiments for HSs are therefore unlikely to be affected by slip. Ye t, with many recen t 



studies of nonlinear c olloidal flow focusing on microscopic properties 



Cohen et al. 



Besseling et al. 



J2006I); 



(120071 )] via microscopy on RI matched suspensions, slip is an important in- 
gredient, and we have shown here that the slip response can be related semi- quantitatively 
to a bulk property of the suspension (osmotic pressure). For other yield stress fluids, in par- 
ticular jammed emuls ions, the slip behavio r may follow a diffe rent powerlaw, i.e. m ~ 0.5 



has been observed in 



Meeker et al. 



(12004H ): 



Seth et al. 



(120081 )] due to elastohydrodynamic 



lubrication for non-repulsive smooth walls. With repulsion, m ~ 1 is recovered, implying 
that our model of the slip-yield transition may carry over to emuls ions. Furthe r ; emu lsions 



with 4> g < 



< 



'rep 



exhibit HS like (Brownian) glassy behavior Gang et al\ (119991 )]. for 



which we thus expect similar Bingham slip behavior as for the HSs. 
In non-Brownian suspensions 



Jana et al. 



(119951 ): ISoltani and Yilmazeri (120081 ): 



Kalyon 



( 120051 )]. slip is also characte rized by m ^ 1 , i.e. a oc v x , i.e. a l ubrica tion layer with a 



thickness independent of v s 



Kalyonl ( 120051 ): lYilmazer and Kalyonl ( 119891 )]. Here, contrary 



to colloids, no slip stress is observed, as expected from the (near) absence of osmotic or wall- 
interaction effects. However, the detailed mechanism for slip in non-Brownian systems is still 
unclear as shown by the different ph enomenolog i cal re 



trated but Newtonian suspensions 



Jana et al. 



ations found: 5/a ~ 0.125 for concen- 
19951)1 .^-inde pendent) , 5/a ~ 0.06 - 0.15 



or pastes of po lydisperse spheres [Soltani and Yilmazeri (120081 )]). and 5/a ~ 2 / [1 — (<p / 4> rcp )] 



Kalyonl (120051 )]. for systems including polydisperse and non-spherical particles). Interest- 



ingly, the latter is similar to our Eq. [IT] for the slip of colloidal glasses, but lacks a theoretical 
basis for non-Bronwian systems. Moreover, in the latter two cases, slip was measured for 
non-uniform stress, such that shear induced migration may affect the interpretation. 

We can also com pare the results with those for (depletion) flocculated colloidal gels in 



Buscall et al\ ( 1l993l )]. There, a linear slip response has also been observed, with 5 decreasing 
from ~ lyum to ~ 10 nm from <fi ~ 0.2 to </> ~ 0.55, without siginificant dependence 
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on particle size or colloid attraction strength. Here the nature of the slip layer is likely 
determined by the </> dependent aggregate lengthscale and structure, rather than the particle 
size. It is worth noting again that in these systems, even significant colloid-wall attraction 
is generally unable to suppress slip; unless the wall roughness is very large, the boundary 
typically acts as a weak 'fracture' plane and a slip response is induced. Insight in the nature 
of this behavior and a theoretical understanding are still lacking. 

In conclusion, we have shown the existence of Bingham-type slip response in colloidal 
HS glasses near smooth non-stick walls. A phenomenological model quantitatively accounts 
for the global rheology and local flow profiles. Slip in HS glasses is effectively caused by 
Brownian motion, creating a lubrication layer and slip response governed by the suspensions 
osmotic pressure, evidenced by the particle size dependence and divergence for — > (p rcp of 
the slip stress and slip viscosity. For HSs, slip is suppressed by colloidal scale wall roughness 
or sufficient vdW wall attraction. Slip can also occur in concentrated liquid-like suspensions, 
but is partial (j a > 7 ^ 0) due to absence of a yield stress. This is also described within 
the phenomenological model, but has only li mited effect on th e bulk rheology. Our study of 



HSs and the study of [Meeker et al\ (j2004bl ); ISeth et al\ (120081 )] for emulsions, together with 



future similar studies for other yield stress fluids, should provide improved predictability of 
yield stress fluid flows in industrial processing and applications. 
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Appendix A: Local and global rheology in cone-plate geometry 



1. Slip at one surface 

This is the case in most of our experiments. Slip occurs at the bottom plate and we set 
o\ = oo (i.e. V\ = 0) and a 2 < a y . In a cone-plate, the relative stress inhomogeneity over 
the gap is ~ 9 2 . This is negligible in our case and we take a uniform at a given r. Using 
Eqs. ( )2P|) with h = rtan(#) we have: 

a(r) =/3rtan(#)[7 a -7(r)] + <x 2 . (Al) 

The stress measured with the rheometer follows from Eq. (Q in Sec. IIVI We define the 
critical applied shear rate 7 ajC = (a y — a 2 )j f3r c t&n(9) such that for 7 a < 7 ac , the bulk shear 
rate vanishes over the entire geometry (regime I). The measured stress is then given by 
Eq. (JTUD in Sec. HVBl When 

7a > 7a, c (regime II), we define r y — (a y — a 2 )/(f3t&n(8) A f a ) such 
that 7 = for r < r y and 7 > for r > r y . The integral in Eq. f fTUj) consists of two parts: 
solid body rotation for r < r y and slip and shear for r > r y . Using Eq. ([5]) and Eq. ([6]) with 
h = rtan(^) and Eq. ( lAlj) gives the following result: 



= ( A2 ) 



1 fl / 3A£ \ ( I _ Aa \ 
T fffi Ur s T 2 




with erg = /3 tan(6 l )r c 7 a , an = aVTa", and Act = a y — a 2 . 

We have assumed that in regime II shear occurs only in the z direction, i.e d r v = 0. But 
because the (partial) slip velocity depends on r, in general d r v 7^ during shear flow. This 
radial velocity gradient can be calculated from v = 'yz, where 7 is given by Eq. ([5]) with 
h = rtan(#), with the result: 



d r v = zdr'j = - 



Act . 2 70 + (2 7a - ^ 



\ \l 1 + i{^ a ~wj / 



(A3) 



Hence, d r v tends to for r — > 00 and r — > r y . However, approaching the boundary of the 
region with solid body, r — > r y , the relative contribution (<9 r t>)/7 grows as: 

dr.v/'j ~ 2z/(r - r y ) < 2rtan(6 l )/(r - r y ). (A4) 
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The shear rate in the vorticity direction may thus be important for r ~ r y , but is negligible 
for (r/r y ) — l ^> 2 tan(#). For our case (9 = 1°), one can neglect d r v in practice, as confirmed 
by the agreement with the experiments, e.g. in Fig. EJ 

2. Slip at both surfaces 

We now also allow slip at the cone with cr 2 < a± < a y . We first discuss solid body 
rotation in regime lb, i.e. o\ < o < a y . The suspensions angular velocity CoVik in this 
regime is determined as follows. The slip velocities at the top and bottom plates are given by 
v\ = (cj— Wbuik)^ and v 2 = Wbuik?", respectively, with u — j a tan(6 l ). Since the total, integrated 
stress a m at the top (a top = ax + 2/3 (u; - Wbuik)r c /3) and bottom {<J bottom = a 2 + 2(3u hn i k r c /3) 
surfaces are equal, the solid body rotation velocity in regime lb is: 

v M k = r7 a tan(fl)/2 + 3(r/r c )(<r 1 - a 2 )/4/3. (A5) 

As a result, for slip at both surfaces, the local stress a(r) differs between the cone and the 
plate, by an amount Aa(r) = (cti — a 2 ) ^1 — JrJ . The rate 7* where the transition from 
slip at the plate (I a ) to slip at the plate and cone (lb) occurs, is determined by equating the 
measured stress due to slip at the plate and that for slip at both surfaces, Eqs. fllQ|l lj) : 

7: = 3(0-i-a 2 )/2/3r c tan(0). (A6) 

When 7 a increases and regime II is approached, the above calculation shows that shear will 
first occur at the edge of the bottom plate for 7 ae = 2%-^un(ff) • ^ ne l° ca l st ress difference 
between the top and bottom plate makes a complete description of shear propagation into the 
cell very difficult. However since a\ — a 2 is often small compared to a y some approximations 
can be made. We assume that shear occurs at the same position r y for the cone and the plate 
and that in the sheared region of the sample (r > r y ) the stress a(z) is uniform. For infinite 
parallel plates we already found that the flow with slip at both surfaces is described in the 
same way as slip at one surface by replacing j5 —>■/?/ 2 and a 2 — > (a% + a 2 )/2, see Eqs. ( 15117]) . 
The same change of variables applied to the cone-plate leads to r y = {2a y — a\ — cr 2 )/f3uj 
and a transition from regime lb to II at 7 aiC = 2 ^~^^ 2 , which, with the approximation 
o\ — o 2 ^ <j y , is equivalent to 7 a>e above. We then again have shear and slip for r > r y and 
solid body rotation for r < r y . The velocity for r < r y follows, as before, from the balance 
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between total stress on the bottom plate and cone: 



rj a tan(0) 



J r<r v 



1 + 



01 - (T 2 



(A7) 



2 \ 2 2<j y — G\ — &2 

For % > 7 a>c , the measured stress is still given by eq. ( 1A2|) with the following changes: 
/3 /3/2 and a 2 {a x + a 2 )/2. 



3. Slip length 

The relation between bulk and applied shear rate in regime II can also be expressed via 
a slip length. For slip at the bottom plate only, we have: 



Is = Vs/^f = /i(7a-7)/7 



(A8) 



with v s the slip velocity, h the gap and 7 and 7 a the local and applied shear rate, respectively 
The full dependence l s (j a ) follows directly from the relation 7(7,2), Eq- © and Eq. §6§ in 



Sect. HVi For large shear rate, Taylor expansion of Eq. §5$) in j a 1//2 results in 



/u/^+oa- 1 ). 

7a 



(A9) 



Figure [19] shows both the exact form for l s from Eqs. ( 1A8II51) and the high shear approxima- 
tion, showing the divergence of l s for % I 7 y and the asymptotic decrease l s ~ 1/^/7^. 

Appendix B: van der Waals interactions 



The van der Waals interact ions between tw o particle s and b etween a particle and a wall 



were calculated according to 



Gregoryl ( I198ll ); iHunterl (l200lh ] . For a particle of radius a 



separated by a distance 5 from the wall we have: 



V, 



pw 



1 + 



5 



5, 
+ - In 



8 



65 V 2a + 5 a \2a + 5. 
while the interaction between two colloids separated by a distance 25 is: 



(Bl) 



V, 



Ai 3 ia 



pp 



2a 



+ 



2a5 



45, (2a + 5 
+ — In ( 5- 

a 



(B2) 



245 \2a + 5 ' (a + 5) 2 a \ (a + 5) 2 , 

For the interaction between the glass plate and particles, the Hamaker constant A 132 , 
where the indexes i = 1,2,3 refer to PMMA, glass, and solvent, respectively (table HT|) . is 
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FIG. 19. Slip length l s versus applied rate for h = 50 /im, a = 10 Pa-s 1//2 , j3 = 1.25 • 10 5 Pa-s-m -1 , 
o y = 10.6 Pa and a s = 2.39 Pa. Full line is the exact form using Eq. (IA8I) . Eq. (J5]) and Eq. ([6]) in 
Sect. IIVI dotted line represents Eq. 



approximately Lee and Komarnenil (120051 )]: 

?>hv e 



= ^ n 



e 3 



+ 



4 fiUW 8v / 2ELiV / ^f + 
while the Hamaker constant between the steel cone (denoted 



3 i=l 



nj — n% 



(B3) 



nr. 



[1) through solvent (3) is estimated using the approach of [Lipkin et all ( 119971 )]: 



by subscript 



2) and the particles 



A 



l:-!2 



3k B T ei - e 3 
4 ei + e 3 



3hu e ni(ni — n 3 ) 



2ni ,/ 2 K + 2 + 2v^2) 



(B4) 



Here is the dielectric permittivity, rij the respective indexes of refraction, h the Planck 
constant, and v e a characteristic frequency which for simplicity is equat ed with the lowest 



energy adsorption peak of PMMA in the ultraviolet v e = 0.3 x 10 16 Hz (jHough and White 
Jl980j)). 
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